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1 Introduction 



Quantum chromodynamics describes the phenomenology of the strong interactions. Currently 
there is a need to improve our knowledge of the higher order structure of the field theory beyond 
earlier one and two loop analysis. Such three loop calculations are performed using efficient 
computer algebra programmes which handle the tedious algebra of such computations. As these 
are computer based it is crucial that independent checks are carried out to confirm the final 
results. We report here on a method which achieves this which is the large Nf self-consistency 
approach developed in Q and applied to four dimensional gauge theories including QED and 
QCD.[^] Essentially the method formulates perturbation theory in an alternative way. For 
models which possess parameters additional to the coupling, such as the number of fundamental 
fields like Nf quarks in QCD, another expansion parameter exists. If Nf is large then 1/Nf is 
small and satisfies the criterion for doing perturbative calculations. Further in this approach 
one considers the critical behaviour of the field theory in the neighbourhood of the d-dimensional 
Gaussian fixed point, where there is a conformal symmetry. As Green's functions obey a simple 
scaling form one computes the associated critical exponent. This encodes information on the 
renormalization of the original Green's function. Then the expansion of the exponent in powers 
of e, d = 4 — 2e, gives the coefficients of the corresponding renormalization group function or 
physical quantity, at that order in 1/Nf. We report here on the application of the technique to 
the renormalization of the physical operators which arise in the operator product expansion of 
deep inelastic scattering. In particular, we focus on the twist 2 flavour non-singlet and singlet 
operators which dominate in most momentum regimes. 3-10 We also discuss the development of 
the method to compute 0{\/Nf) information on the process dependent moments of the structure 
functions based on the earlier work of Broadhurst and Kataev[ll, 12] inspired by El. 



2 Formalism 

The twist 2 Wilson operators whose anomalous dimensions we compute in critical exponent form 
are! 
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D^ n q — trace terms 
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where q is the quark field, I < i < Nf , 1 < I < N c , n is the moment 

f^A^A^, A^ is the gluon field, 1 < a < N c , the S denotes symmetrization of the Lorentz indices 
and Tfj are the group generators. We discuss the technique to deduce the anomalous dimensions 
of (l)-(3) by considering the non-singlet case first. |jl4| To carry out a perturbative analysis one 
inserts the operator into some Green's function and determines the pole structure with respect to 
some regularization. In the critical point approach^ one deduces critical exponents by using, 
instead of the usual propagators, the critical ones whose structure is determined solely from 
scaling arguments and dimensional analysis. H, [| So in d-dimensions the asymptotic scaling 
forms of the quark and gluon propagators in the critical region, k 2 — > oo, are 

Aft 
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where d = 2fi. The fixed point is defined to be g c where f3(g c ) = 0, g c ^ 0. From the one loop 
/3-function in rf-dimensions fl6|, O, 111] 



(i{g) = {d-A)g + 



2 -T(R)N f - ^C 2 (G) 



g 2 +0(g 3 ) (2.6) 



where Tr(r a T b ) = T(R)5 ab , T a T a = C 2 (R)I, f acd f bcd = C 2 (G)5 ab and g = (e/2vr) 2 is our 
dimensionless coupling constant. Then from (6) at leading order, as Nf — > oo, 

g c = — ■ ()\ - } (2.7) 




T{R)N f 

The dimension of the fields in (4) and (5) are determined from the fact that the action is 
dimensionless and their anomalous pieces are defined via 

a = M-l + fr, , = 1-rj-x (2.8) 

where rj is the quark anomalous dimension and x is the anomalous dimension of the quark gluon 
vertex and each have been calculated at 0(1/Nf) in the Landau gauge, [g] The quantities A and 
B in (4) and (5) are the amplitudes of the fields and b is the covariant gauge parameter. As 
(l)-(3) are physical then their anomalous dimensions are gauge independent. || So calculating 
in an arbitrary gauge and observing the cancellation of b provides a non-trivial check on any 
calculation, though the work of || Q was carried out in the Feynman gauge. 

To proceed in the critical point analysis the regularization is introduced by shifting /?—>/? 
— A.JHJ Then the residue of the pole in A when (1) is inserted in a Green's function (qO^q) 
contributes to the renormalization or exponent of the bare operator, 7£^(g c )- The full gauge 
independent exponent is given by 

7^(5c) = r, + 7^ (9c) (2.9) 
Then we find the anomalous dimension exponent of (1) is 

(n) 2C 2 (R)<ja - l) 2 V? \ (n-l)(2/i + n-2) 



(2n-l)(p-2)T(R) 
+ -^rr^-l + n)-^)] 



(ji + n-l)(p + n-2) 



(2.10) 



in arbitrary dimensions, where Tp{x) is the logarithmic derivative of the r-function and the 
subscript 1 denotes the coefficient oil/Nf in the expansion of 7^g(5c)- 

There are several checks on (10). First, when n = 1, the original operator corresponds to 
a symmetry generator and therefore its anomalous dimension ought to vanish [ 20 1. It is trivial 
to see that 7^g i(gc) = 0. Second, performing the e-expansion of (10) and comparing with the 

explicit perturbative function 7jjg(<7), evaluated at (7), the coefficients agree with the leading 
order two loop analytic forms given in[]|, |], Recently information on the full 3-loop structure 
has been provided. pl| That calculation made extensive use of a computer algebra programme 
to the extent that 7ns(50 ^ s known at 3-loops for the even moments up to n = 10. Expanding 
(10) to 0(e 3 ) one obtains the analytic expression 

= ls,{n)-f 7 S 2 {n)-^S 1 {n) + f 2 
[12n 4 + 2n 3 - I2n 2 - 2n + 3] 

(2.11) 



27n 3 (n + l) 3 

for the 3-loop coefficient at 0(1/Nf), where Si(n) = Ya=i Substituting for the various n 

we record that (11) is in exact agreement with^TJ. Analytic expressions at 0(1/ Nf) can be 
obtained for the higher order coefficients. E3] 
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3 Singlet Anomalous Dimensions 



We now turn to the singlet case. As (2) and (3) have the same dimensions they mix under 
renormalization which complicates the computation of their anomalous dimensions. By con- 
trast to (1) one has a matrix of anomalous dimensions, 7^ (<7)-|| Its eigenvalues correspond 
to the anomalous dimensions of two independent combinations of the original (bare) operators 
which are of physical interest. Perturbatively 7«^(<?) is determined in the same way as 7^3(5) 
but the operators are also included in a gluon 2-point function (AOA) and it, and therefore 
its eigenvalues -y±\g), have been computed to two loops. 7-10 To construct the corresponding 
anomalous dimension exponents, one analyses the leading order large Nf graphs at criticality as 
discussed before but now determines the matrix of critical exponents. The graphs for these have 
been given elsewhere |7| but we note that there are several two loop graphs which contribute at 
leading order in l/Nf. We find the eigenvalues of 7}" (g) at criticality for n even, at leading 
order, are 

7 i n) (<7c) = -2(/x-2) (3.1) 

_ 2(^-l) 2 (n-l)(2 / Li + n-2) 



(fi + n- l)(n + n- 2) 



(2n-l)0M-2)T(R)N f 
+ 4m(m-1)[^(a*-1 + »*)-^(m)] 

- i)r(n - i)r(2/i) 

(ji + n-l)(ji + n- 2)r(2/x - 1 + n) 
x [(n(n - 1) + 2(/i - 1 + n)f 

+ 2{jjl - 2)(n(n - l)(2/i - 3 + 2n) + 2{fi - 1 + n))) 



(3.2) 



In the critical point approach (12) and (13) follow naturally without diagonalization as 7^ (g c ) 
is triangular at leading order. Eq. (12) represents the anomalous dimension of the operator 
which has (3) as its dominant contribution in the linear combination we mentioned, whilst (13) 
corresponds to (2) being dominant. The Nf dependence of each term is different due to the 
particular way in which quark loops occur at leading order in the Green's functions. Specifically 
the e-expansion of (12) involves only a one loop term and no subsequent terms, consistent with 
the perturbative analysis. (We also note that recently the renormalization of the gluonic sector 
has been re-examined to understand the role infrared divergences play.p2[) 



There are several checks on (13). When n = 2, the operator corresponds to the energy mo- 
mentum tensor and since it is conserved its anomalous dimension vanishes. [p0|] Clearly, j_ '(c/ c ) 
= 0. Secondly, using the explicit 2-loop matrix, 7}" (g), and extracting both eigenvalues as a 
perturbative expansion in g, the coefficients of (13) to 0(e 2 ) agree exactly. 7-10 We have also 
checked the cancellation of the gauge parameter. For comparison sake we record that the three 
loop leading order 1/Nf coefficient of 7_ (ff) is 

V> 10 ,,u 17 2(n 2 + n + 2) 2 [S 2 (n) + S 2 1 (n)} 

«3 = gS 3 (n) - -S 2 (n) + - - - 2)(n - 1)2(n _ 1} 

- 2S 1 (n)[n 9 + 6n 8 - 36n 7 - 216n 6 - 552n 5 - 810n 4 - 81 In 3 
- 690n 2 - 132n + 72]/[27(n + 2) 2 (n + l) 3 (n - l)n 3 ] 

- [100n 10 + 682n 9 + 2079n 8 + 3377n 7 + 3389n 6 + 3545n 5 + 3130n 4 

+ 118n 3 - 940n 2 - 72n + 144]/[27(n + 2) 3 (n + l) 4 n 4 (n - 1)] (3.3) 
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4 Structure Function Moments 



A second ingredient in the Wilson expansion is the determination of the process dependent 
Wilson coefficients, C l (q 2 /m 2 , g), i = NS, ip or G, whose momentum evolution is controlled by 
the anomalous dimensions. As these have also been computed at 0(g s ) and to n = 10 in pi] 
it is equally important to have an efficient method of calculation in large Nf. By contrast with 
(10,12,13), which relate to the divergent part of Green's functions, the C l 's are determined from 
the finite part of the amplitude after renormalization. The most efficient way to achieve this is 
based on [11, [l^] where one computes a minimal set of one or two loop graphs. In that approach 



the relevant diagram is calculated in strictly four dimensions but where the gluon line has an 
exponent of 1 + 5 instead of 1. Multiplying the resultant expression, which will be a function 
of 5, by e 55//3 the coefficients of its Taylor series in 5 correspond to the perturbative coefficients 
of C % at 0(1/Nf). This exponential factor is necessary to ensure results are in the MS scheme, 
and the 5/3 arises from the finite part of the charge renormalization. (We note that we have 



developed this independently of a similar approach in [23] for a different problem.) We have 
applied this method to the determination of the simplest process ie the non-singlet longitudinal 
Wilson coefficient in the Bjorken limit. Including the crossed amplitude we find, at the Lth 
loop, 



NS 
long 



d5 L 



8C 2 {R)e 5S / 3 r{n + 5)g 



(2 -S)(l- 5){n + 1 + <5)r(n)r(l + 5)x n 



(4.1) 



5 = ±N f T(R)g 



We recall that the coupling constant expansion of C£^(l,5) is O(g) and not 0(1). Expanding 
in powers of 5 the coefficients agree with the 2-loop results of [Q, ^4|, |25[] for all n. We note that 
the three loop coefficient at 0(1/Nf) is 



NS 
"long,3 



16C 2 (R)T 2 (R) 
9(n + 1) 

+ 2Si(rO(H 



^ - S 2 (n) + S 2 (n) 
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n 



1 



n+1 



+ 



+ 



19(2n + 1) 
3n(n + 1) 
2 

+ 



n- 



(n + 1) 2 n(n+l) 



(4.2) 



which agrees with the even moments of [21] up to n = 10. 



5 Conclusions 



We have presented the latest results of applying the large Nf self-consistency programme to the 
operator product expansion of deep inelastic scattering. Analytic results have been obtained 
for the higher order structure of the anomalous dimensions of the physical operators and the 
moments of the structure functions at 0(1/Nf). An important feature of our results is that we 
have given an insight into the (complicated) analytic structure at three and higher loops as a 
function of n, primarily by exploiting the conformal symmetry of the d-dimensional fixed point. 
This is a first step in gaining an insight into the subleading 0(1/ N 2 ) coefficients in relation to 
determining the 0(g 3 ) Altarelli-Parisi splitting functions. Further, we have indicated how to 
deduce useful information compactly and efficiently, for the moments of Wilson coefficients in 
1/Nf. Although we have focussed on the simplest case, the longitudinal part of the non-singlet 
amplitude, it ought to be possible to adapt the approach to determine information on the singlet 
structure. 
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